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ABSTRACT-Whereas Cantor multilayers made of an isotropic dielectric-magnetic mate- 
rial with positive refractive index will show power-law characteristics, low-order Cantor 
multilayers made of materials with negative refractive index will not exhibit the power- 
law nature. A reason for this anomalous behavior is presented. 
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1 Introduction 

This letter addresses the incorporation of isotropic materials with negative refractive 
index M, 0] in fractal filters inspired by Cantor dusts H, BJ. 

The emergence of Cantor dusts, bars and cakes during the late 19th century has 

been described at some length by Mandelbrot ||. Briefly, the simplest Cantor dust 

is formed by dividing the closed interval [0, 1] into 3 pieces and removing the center 

open piece (1/3, 2/3), repeating the trifurcation-and-removal process on the remaining 

intervals [0, 1/3] and [2/3, 1], and continuing in that fashion ad infinitum. The fractal 
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(similarity) dimension of the resulting dust is log 2/ log 3 ~ 0.6309. Similar structures 
in p-dimensional space, (p = 1, 2, ...), can be constructed via spatial convolution Q. In 
particular, the constructs called Cantor bars appear to have captured the imagination of 
optical-filter researchers, as recounted recently by Lehman ||, because of their putatively 
self-similar response properties in the frequency domain Q. 

The materials of choice for optical Cantor filters are isotropic dielectric with relative 
permittivity e r . Although e r is a complex-valued function of frequency, the usual practice 
in optics is to ignore dissipation by setting Im [e r ] = 0. In the area of fractal optics, with 
emphasis still on understanding basic interactions in nonperiodic multilayers, dispersion 
is also ignored ||— (§[]■ The structural self-similarity of the Cantor bars is then expected to 
result in the self-similarity of the spectral reflectance/transmittance responses of optical 
Cantor filters to normally incident light §. Truly, physically realizable Cantor filters 
are not actually fractal but pre-fractal instead || - - so that the spectral self-similarity 
can only be approximate |l( 



On examining the available literature, two questions arise. First, will the situation 
change for Cantor filters made of isotropic dielectric-magnetic materials (with relative 
permeability denoted by \i r > 1)? Second, will the situation change if both e r < and 

fir < 0? 

The second question arose because of the supposed verification of the existence of 
negative refractive index (NRI) by Shelby et al. last year. Experiments performed on 
certain composite materials with oriented microstructure suggested that these materials 
are endowed with negligible dissipation as well as NRI in some appreciably wide fre- 
quency band in the centimeter-wave regime. Also called left-handed materials by some 
researchers (despite possessing no handedness), in NRI materials the phase velocity is 
pointed opposite to the direction of energy flow (and attenuation) Q, jy]]. Although 
the extant experimental results are not perfect fll2| ) jDJ, the essential conclusion of the 
existence of NRIs appears undeniable. As NRIs can potentially lead to exciting new 



technologies [14], theoretical consideration is warranted. 
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In this letter, we answer the two questions posed earlier in a unified way. Section 2 
is devoted to the theory of reflection and transmission of normally incident plane waves 
by Cantor multilayers. Numerical results are presented and discussed in Section 3. 



2 Theory 

A Cantor multilayer is constructed sequentially as follows: Take a layer of thickness £q 
made of a certain material with e r and p r as its constitutive parameters. Call this layer 
a multilayer of order N = 0. Next, cascade two multilayers of order N = inserting a 
space of thickness £q/ f, f > 1, in between. Call this a multilayer of order N = 1. Its total 
thickness £\ = (2 + 1//)£q. Continue in this manner. Thus, a multilayer of order N + 1 is 
formed by inserting a space of thickness £n / f between two multilayers of order N. The 
thickness of a multilayer of order A + 1 is then £ N+1 = (2 + l/f)£ N = (2 + l/f) N+1 e . 
The fractal dimension of the multilayer is given by 

V= , (1) 

log(2 + l//)' U 

which concept is applicable strictly in the limit iV — > oo. 

Let a Cantor multilayer of order N occupy the space < z < £jy. Suppose a plane 
wave is normally incident on this multilayer from the vacuous half-space z < 0, with Ao 
denoting its wavelength. Therefore, a reflected plane wave also exists in the same half- 
space. Furthermore, a transmitted plane wave is engendered in the vacuous half-space 
z> In- The corresponding electric field phasors are given by 



E(z) = < 



exp(ifco-z) + Pn exp(— ik^z) , z < 

Tjvexp [iko(z — £n)] , z > £n 



where ko = 2tt/\q is the wavenumber in vacuum; and tn are the reflection and 
the transmission coefficients, respectively, both complex-valued; and (w^,^,^) is the 
triad of cartesian unit vectors. An exp(— iujt) time-dependence is implicit, where to = 
ko/^eopo) 1 / 2 is the angular frequency, while eo and po are the permittivity and the 
permeability of vacuum, respectively. 



The coefficients p^ and tat can be easily determined using a 2 x 2 matrix algebra 
After defining the two matrixes 



A 



the matrixes M p , < p < N, are iteratively computed as 



M p+ i = M p • e ^ p//)A • M p , < p < N - 1 






Mo 


B = 
















eoe r 



(3) 



beginning with 



1 = e 



iu)£o B 



(4) 



(5) 



The boundary value problem for the electromagnetic fields then involves the solution of 
the equation 





1 


= Miv • 



















0- + PN) 

V^i 1 ~ Pn) 



(6) 



where ijq = (/io/eo) 1//2 is the intrinsic impedance of vacuum. The principle of conservation 
of energy entails that |pat| 2 + |tat| 2 < 1, with the equality coming in when the multilayer 
is made of a non-dissipative material. 



3 Numerical results and discussion 



Following normal practice, we implemented the foregoing equations to compute pn and 
ttv for non-dissipative and non-dispersive materials. We varied the quantity £ = fco^o 
for various values of N, while keeping e r and fi r fixed. 

The spectrums of \pn\ 2 and |ttv| 2 turned to be identical to the ones reported in the 
literature 0, || for optical Cantor filters (i.e., with e r > 1 and p r = 1). Those for 
{e r > 1, p, r > 1} and {e r < 0, fi r < 0} turned to be qualitatively similar, and therefore 
do not need reproduction here. 

As £ increases from zero, the fundamental layer thickness £q becomes an increas- 
ingly significant fraction of the wavelength Aq = 2Tr/ko, and eventually surpasses Aq. In 
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other words, layers are electrically thin for small £, and an increase in £ amounts to 
magnification. Therefore we evaluated the value (jv of C at which the first minimum 
of |tjv| occurs as £ increases from zero, thereby reckoning (jy as a reasonable parameter 
containing structural information on the chosen multilayers. If indeed the structural self- 
similarity of Cantor multilayers would result in their spectral self-similarity, we expect 
the relationship 

Gv = 2-^Co (7) 
to emerge from our numerical investigations. 

Figures 1 and 2 contain plots of log (jv versus N for Cantor multilayers made with 
positive refractive index (PRI) materials ({e r = 3, \x T = 1.02} or {e r = 4, fx r = 1.02}), 
and for Cantor multilayers made with their NRI analogs ({e r = —3, \x r = —1.02} or 
{e r = —4, n r = —1.02}). The factor / = 1 for Figure 1, and / = 2 for Figure 2. 
Two conclusions can be immediately drawn from these two figures as follows: 

A. The relationship (jv = 2~ N l' Dpm satisfied by Cantor multilayers with PRI ma- 



terials is a power law with Vpm > T>, and could be fractalesque |TH, 17 



B. The data for Cantor multilayers with NRI materials indicates two different regimes, 
one for small N and the other for large N, the second regime characterized by a 
power law. 

The foregoing conclusions suggest that the effect of NRI materials on electromagnetic 
fields must be substantively different from that of PRI materials, for the anomalous first 
regime to arise for Cantor multilayers with NRI materials. Furthermore, in the present 
context, the difference must be evident definitely for order iV = 0. 

Hence, we analyzed the planewave response of a single layer to obtain 

= (^r ~ l)sin/? 

Po (»^ + l)sin/3 + 2i^cos/3 1 ' 



and 



T ° (»/2 + l)sin/3 + 2if7 r cos/3 ' {> 



Here, the relative impedance rj r = +y p r /e r must be positive real, while the sign of 
j3 = ko£Q^/fi r e r has to be positive/negative for PRI/NRI materials p^| . Denoting the 
phase of a complex number £ by Z£, we conclude from the foregoing equations that 

{e r -e r ,/JL r -> -p r } =>• {\p \ -> \po\, \t \ -> |r |, Zp -> -Zp ,/T -> -Zr } . (10) 



In light of the relationship (10), let us compare a PRI layer and a NRI layer - 
labeled a and 6, respectively — such that p Ta = —p rb > and e Ta = —e rb > 0, while the 
wavenumber ko is fixed. If the thicknesses of the two layers are such that the sum /? a + |A>| 
is an integral multiple of 2ir, then (|8j) and @ yield po a = po b and ro Q = r 0b . Thus, a PRI 
layer of a certain thickness is equivalent to a NRI layer of different thickness, in terms of 
the complex-valued reflection and transmission coefficients at a fixed wavelength. But 
the thickness of the equivalent NRI layer is wavelength-dependent — which implies that 
a PRI Cantor multilayer is equivalent at different wavelengths to different NRI Cantor 
multilayers. Not surprisingly therefore, the spectral characteristics of a PRI and a NRI 
Cantor multilayers with the same l§ are not isomorphic. 

The difference is very noticeable for small N in Figures 1 and 2. As N increases, 
the value of Cat decreases for both PRI and NRI multilayers — in other words, the 
fundamental layer of thickness 1$ becomes electrically thinner at the first transmittance 
minimum and, therefore, a weaker reflector as well as a stronger transmitter. Structural 
characteristics then dominate over the consequences of (|10"D, because \po\ ~ and |ro| — 1. 
As the difference between PRI and NRI multilayers lessens with increasing N, the latter 
also begin to evince power-law characteristics. 

The crossover between the anomalous and the power-law regimes for NRI multilayers 
takes place at a higher value of N as / increases. This general trend is indicated by 
Figures 1 and 2 as well as calculations for other values of /. 

To conclude, we have shown that the planewave reflection and transmission spec- 
trums of a Cantor multilayer made of an isotropic dielectric-magnetic material with 
positive refractive index shows power-law characteristics which indicate spectral self- 
similarity. However, if the same multilayer were to be made of a material with negative 
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refractive index, then the power-law nature is not going to be evident when the interac- 
tion between the material layers and the interleaving vacuous spaces is substantial (the 
small- N regime). The existence of this anomalous regime can be attributed to the re- 
flection/transmission phase reversal of a NRI layer in relation to its PRI analog. If that 
interaction is insubstantial (the large-iV regime) , the structural features would dominate 
the constitutive features, and the power-law characteristics would be evident also for the 
NRI Cantor multilayer. 

Acknowledgement We thank an anonymous reviewer for the suggestion to examine 
the large-iV responses. 
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Figure 1: Calculated values of log Cat for Cantor multilayers of orders N when / = 1. 
Dotted lines join points for PRI Cantor multilayers, dashed-dotted lines for NRI Cantor 
multilayers, and solid lines for (jv = 2~ n ' t> ^q. (a) e r = ±3 and fi r = ±1.02; (b) e r = ±4 
and fj, r = ±1.02. 
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Figure 2: Same as Figure 1, but for / = 2. 



